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ABSTRACT 

We compare the cosmological first-order post-Newtonian (1PN) approximation with the relativistic 
cosmological linear perturbation theory in a zero-pressure medium with the cosmological constant. 
We compare equations and solutions in several different gauge conditions available in both methods. 
In the PN method we have perturbation equations for density, velocity and gravitational potential 
independently of the gauge condition to IPN order. However, correspondences with these IPN equa- 
tions are available only in certain gauge conditions in the perturbation theory. Equations of perturbed 
velocity and the perturbed gravitational potential in the zero-shear gauge exactly coincide with the 
Newtonian equations which remain valid even to IPN order (the same is true for perturbed velocity 
in the comoving gauge) , and equations of perturbed density in the zero-shear gauge and the uniform- 
expansion gauge coincide to IPN order. We identify other correspondences available in different gauge 
conditions of the perturbation theory. 

Subject headings: cosmology: theory — large scale structure of universe 



1. INTRODUCTION 

Einstein's gravity is generally accepted as the gravity 
to handle astronomical phenomena. The theory holds 
a remarkable track record in the solar-system test based 
on vacuum Schwarzschild solution and the parameterized 
post-Newtonian approximation where the gravitational 
fields are supposed to be weak. It is true that Einstein's 
theory has not failed in any experimental test based on 
modern scientific and technological development up till 
today, but it is also true that there has been no experi- 
mental test of the theory in the strong gravitational field 
and in large scale even including the galactic scale. Ein- 
stein's gravity is generally accepted in cosmology mainly 
based on its successes in other astronomical and Earth 
bound tests and the theory's own prestige associated 
with Einstein's fame and historical legacy. 

A self-consistent treatment of cosmological world 
model is possible in Einstein's gravity. Without a lead by 
Einstein's gravity, however, the spatially homogeneous 
and isotropic cosmological world model based on New- 
ton's gravity is known to be incomplete and indetermi- 
nate (Layzer 1954; Lemons 1988). Despite such troubles 
in the background world model, evolution of perturba- 
tions in Newton's theory is known to be quite success- 
ful in reproducing the corresponding results in Einstein's 
gravity (Lifshitz 1946; Bonnor 1957; Noh & Hwang 2004; 
Jeong et al. 2011). Considering the action- at- a-distance 
nature of Newton's gravity, such a coincidence is a non- 
trivial result. Differences between the two theories, how- 
ever, appear as the scale approaches the horizon. We will 
address this issue in this work. 

If we accept Einstein's theory in analyzing the large- 
scale cosmic structure in current era, we have two meth- 
ods available. One well known method is the pertur- 
bation theory where all dimensionless deviations in the 
metric and the energy-momentum tensor from the back- 
ground world model are assumed to be small. If we ac- 
cept only linear order deviations we have the linear per- 



turbation theory. The perturbation theory assumes small 
deviation but is fully relativistic, generally applicable in 
all scales including the super-horizon scale and to par- 
ticles with relativistic velocities (Lifshitz 1946; Harrison 
1967; Nariai 1969; Bardeen 1980; Peebles 1980; Kodama 
& Sasaki 1984; Bardeen 1988; Mukhanov et al. 1992; Ma 
& Bertschinger 1995). 

The other less known method is the cosmological post- 
Newtonian (PN) approximation where all dimensionless 
deviations are assumed to be weakly relativistic with 
<C 1 (M and R are characteristic mass and length 

scales) and for a virialized system % <C 1 (v is the char- 
acteristic velocity involved). The first-order PN (IPN) 

approximation makes expansion up to ~ \ order. 
The PN approximation assumes small relativistic effects 
and is applicable only in sub-horizon scale. But the equa- 
tions derived are applicable to fully nonlinear situations 
(Chandrasekhar 1965; Futamase 1988; Tomita 1988; Fu- 
tamase 1989; Tomita 1991; Futamase 1993a, 1993b; Shi- 
bata & Asada 1995; Asada & Futamase 1997; Tanaka & 
Futamase 1999; Hwang et al. 2008, PN2008 hereafter). 
Therefore, the two methods are complementary to each 
other. 

If we encounter cosmological situations where both 
nonlinearity as well as relativistic effects are important 
we may need full-blown numerical relativity implemented 
in cosmology. Currently such a general relativistic nu- 
merical simulation in cosmology is not available. The 
nonlinear perturbation analysis, being based on pertur- 
bative approach, is not sufficient to handle the gen- 
uine nonlinear aspects of structure formation accompa- 
nied with self-organization and spontaneous formation of 
structures. In order to handle the relativistic nonlinear 
process in cosmology we believe the post-Newtonian ap- 
proach is currently practically relevant to implement in 
numerical simulation. 

We can find cosmological situations where the cosmo- 
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logical post-Newtonian approach, being weakly relativis- 
tic but fully nonlinear, might have important applica- 
tions. Especially, the current cosmological paradigm fa- 
vors a model where the large-scale structures (requir- 
ing the relativistic treatment) are in the linear stage, 
whereas small-scale structures are apparently in fully 
nonlinear stage. The often adopted strategy is to as- 
sume the small-scale nonlinear structures are fully under 
control by the Newtonian gravity. In the galactic and 
cluster scales we have the general relativistic measure 
~ ~ 10~ 6 — 10~ 4 , thus small but nonvanishing, 
and indeed the 1PN (weakly relativistic) assumption is 
quite sufficiently valid. Thus, we believe the 1PN ap- 
proach would be quite relevant to estimate the general 
relativistic effects in the nonlinear clustering processes of 
the galaxy cluster-scale and the large-scale structures. 

In this work we will compare the two relativistic meth- 
ods in the matter dominated era: the 1PN method vs. 
linear perturbation theory. The 1PN approximation is 
based on previous studies in Chandrasekhar (1965) and 
PN2008, whereas the linear perturbation theory is based 
on previous studies in Bardeen (1988), Hwang (1994), 
and Hwang & Noh (1999). We will compare the equa- 
tions and solutions derived in the two methods. In both 
methods we have gauge degrees of freedom which need 
to be fixed by the gauge conditions. The Newtonian per- 
turbation theory appears as the zeroth-order PN (OPN) 
approximation. Thus, we naturally also have Newtonian 
theory for the comparison. To 1PN order we will show 
that the equations for density, velocity and gravitational 
potential do not depend on the gauge conditions with 
each variables gauge-invariant to 1PN order. In the per- 
turbation theory, however, the perturbation variables for 
density, velocity, potential, curvature, and other kine- 
matic variables (like expansion and shear) do depend on 
the gauge conditions adopted. Our emphasis in this work 
is on the correspondences between the PN variables and 
perturbation theory variables based on different gauge 
conditions. 

For the velocity and gravitational potential the New- 
tonian (OPN) equations are valid even to 1PN order, 
whereas for the density variable we have 1PN correction 
terms: see Equations ([25" ]l -([2"8" ]) . The zero-shear gauge 



will be shown to be distinguished by showing the per- 
turbed velocity and gravitational potential having exact 
correspondences with the OPN (thus 1PN as well) equa- 
tions, and the perturbed density showing correct corre- 
spondence with the 1PN equation. The zero-shear gauge 
can be contrasted with the comoving gauge where the 
perturbed density and velocity have exact Newtonian 
correspondences even to the second order perturbations 
in all scales whereas the gravitational potential vanishes 
to the linear order. Other PN correspondences of the 
perturbation theory are summarized in Section [5j 

2. 1PN APPROXIMATION 

To the 1PN order our metric convention is (Chan- 
drasekhar 1965; Chandrasekhar & Nutku 1969; PN2008) 



ds 2 



1 1/, 

-2U+-r 2U 2 



4$ 



c 2 dt 2 



2aP t cdtdx l + a 2 1 + -^2V ^dx'dx 3 , (1) 



where a(t) is a cosmic scale factor. We consider a flat 
background, thus 7^ = <5y. Then, we have V = U 
(PN2008). By ignoring c~ 2 a 2 Cij part in the g+j we al- 
ready have taken spatial gauge conditions without losing 
any generality: see Section 6 in PN2008. 

2.1. Equations in the zero-pressure case 

We consider a zero-pressure fluid in the presence of cos- 
mological constant. The 1PN energy-momentum tensor 
is presented in Equation (21) of PN2008 (Chandrasekhar 
1965). With vanishing pressure, internal energy, flux, 
and anisotropic stress, we have 



T ab = im a u b , [i = qc 2 , 



c a 



(2) 



where \x and u a are the covariant energy density and the 
normalized fluid four-vector, respectively. To the linear 
order, we have (PN2008) 
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Pi), T?=0. (3) 



The 1PN equations are (Chandrasekhar 1965; PN2008) 
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(4) 



(5) 



(6) 



(7) 



where a vertical bar indicates spatial covariant deriva- 
tive based on 7^ as the metric. These follow from the 
mass-conservation, momentum-conservation, Raychaud- 



hury (Go — G-), and momentum-constraint (G°) equa- 
tions, respectively; for the general case, see Equations 
(114), (115), (119) and (120) in PN2008. Terms on the 



POST-NEWTONIAN APPROXIMATION VS. PERTURBATION THEORY 



3 



left-hand-side provide the Newtonian (OPN) limit, and 
the ones in the right-hand-side are 1PN contributions. 
Notice that the 1PN terms include up to fourth-order 
perturbations. 

To perturbed order we have g = Qb + Sg where we will 
ignore the subindex b indicating the background quantity 
unless necessary. 

Our cosmological PN approach assumes near fiat back- 
ground, but is valid in the presence of the cosmological 
constant. Equations for the background 
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0. 



(8) 



were subtracted in deriving the PN equations; see Section 
3.2 in PN2008. 

2.2. Newtonian equations as the OPN limit 

The OPN approximation gives the Newtonian limit. 
The OPN metric is 



ds 2 



1 - -^2U ) c z dt z + cfSijdx'd.r 1 . ( !) ) 



(10) 

(11) 
(12) 



To OPN limit Equations (g}-© give 

S = --V-[(1 + S)v], 
a 

■ a 1 r, l ^rr 

v + -v + - v Vv = - W, 

a a a 

~U = AirGgS, 
a 2 

where 5 = Sg/g, v = v l , and <5$ = —U are the rela- 
tive density contrast, velocity perturbation, and the per- 
turbed Newtonian gravitational potential, respectively. 

In the Newtonian context, assuming the presence of 
spatially homogeneous and isotropic background the 
above equations follow from the mass conservation, the 
momentum conservation, and the Poisson's equation, re- 
spectively (Peebles 1980); for weakness of Newton's the- 
ory in handling background world model, however, see 
Layzer (1954) and Lemons (1988). In fact, the com- 
plete nonlinear Newtonian hydrodynamic equations are 
already built in as the OPN approximation of Einstein's 
gravity (Chandrasekhar 1965; Bertschinger & Hamil- 
ton 1994; Kofman & Pogosyan 1995; see Section 4.2 in 
PN2008). Up to this point is a review of the cosmological 
1PN approach. 

To the linear order, from Equations (ITUl) - (fl"2l we can 
derive 



2HS - AnGgS 



1 



a 2 H 



a 2 H 2 



= 0, (13) 



(aV • v)" + H (aV • v)' - AirGg (aV • v) = 0, (14) 

(all)" + 2H (all)' - AirGg (aU) = 0, (15) 

where H = a/a. Assuming K = = A, the growing 
mode solutions are 

2 C " A r l x, v 2c2A r 1 n 3 r 

(16) 



where C(x) is an integration constant; the normaliza- 
tion of the constant C will become clear later, see be- 
low Equation (|75|) . Notice that as the Newtonian (OPN) 
equations are not aware of the presence of horizon, New- 
tonian equations are supposed be valid in all scales. The 
Newtonian or the OPN order equations are not affected 
by the gauge transformation. 

2.3. Linear limit of 1PN equations 

By keeping only linear order terms Equations (JT]) 
give 



i (a 3 Sg) + = -^eU, 
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(19) 



(20) 



(21) 



By taking a divergence, Equation (|20|) gives 

AnGgav 1 ^ = A ({j + HU} , 

which follows from the OPN order of Equations (|17p and 
(|19[) . By introducing a combination 
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1 r 
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2$ * + (aPi)' 



(22) 



(23) 
(24) 



From Equations ([T7l) - ()19|) we can derive (Takada & 
Futamase 1999) 



Equations (fl"%]) and P^|) become 

- (avi) ~ -Ua = 0, 
a a 

A „„ 3/-- a ■ a 

-^U + AnGSg = — - [U + 3-U + 2-U 



S + 2HS- AnGgS = -| ( -U + 2-U 
c z \ a a 



(25) 



Since terms in RHS are already 1PN order, using Equa- 
tion (|19l) to the Newtonian order, we have 



6 + 2H5 - AnGgS = 



12nGga 2 
c 2 A 



HS 



(2H + H 2 ^j 



Similarly we have 

(aV ■ v)" +H(aV ■ v)' - AirGg (aV • v) = ( 

(aU)" +2H {aUy -ATrGg{aU) = 0, 

where we used Equation (|21[) to appropriate order. The 
terms in the RHS of Equation (|26[) are linear contribu- 
tions from 1PN correction terms. This apparently shows 



S . 

(26) 

(27) 
(28) 
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that the cosmological PN approximation is valid only far 
inside the horizon with c 2 A >• 12irGga 2 ~ H 2 ~ if; 
as we approach the horizon scale higher order PN cor- 
rections, with ( a c /^ )" ~ (^t)™ order correction factors 
for nPN order, become important. 

It is remarkable that for the perturbed velocity and 
perturbed gravitational potential in Equations (|27|) and 
(|28p we have no PN correction terms to the 1PN order. 
We have not imposed the temporal gauge condition to 
derive Equations (l25l) - ([27l) . and the variables 5, v, U 
and hi are gauge-invariant: this is the subject we discuss 
below. 

In our 1PN metric convention in Equation ([T]) we 
already have taken spatial gauge condition by setting 
9ij = a 2 (l + c~ 2 2V)8ij\ for a thorough examination of 
the gauge issue in the PN approximation see Section 6 
in PN2008. Under the remaining gauge transformation 
x a = x a + £ Q (x e ) with 



e = 1^(2)0 + 1 £(4)0 (2g) 
c C J 



we can set £( 2 )° = without losing any generality: see 
Equation (173) in PN2008. Then, U does not depend on 
the gauge and we have [see Equations (175) and (176) in 



PN2008] 

Pi=Pi-\t W ° ti , $ = *+^ (4)0 . (30) 

Thus, in the PN approach we have freedom to impose the 
temporal gauge (hypersurface) condition on P 1 ^ or 

A combination 2$^ + (aPi)' , thus hi, is gauge- invariant. 
The fluid variables Sg and v in our case of vanishing inter- 
nal energy and the flux are gauge invariant: see Section 
6 in PN2008. 

In PN2008 we have introduced several temporal gauge 
conditions summarized in Table 1. In Table 1, the 
propagation speed indicates the speed of propagation 
of the gravitational potential U in the Equation (| 19|) 
or (|24p . For example, the harmonic gauge condition 
makes the Laplacian operator A in the OPN limit re- 
placed by a d'Alembertian operator □ by the 1PN correc- 
tion terms, thus making the Poisson's equation in OPN 
limit to a wave equation with the propagation speed c 
by the 1PN correction. Similarly, the uniform-expansion 
gauge (and the Chandrasekhar's gauge) leaves the action- 
at-a-distance nature of the Poisson's equation, and the 
transverse-shear gauge makes Equation (fl9|) to be no 
longer a wave equation. Apparently, the propagation 
speed of the gravitational potential depends on the gauge 
choice. However the propagation speed of covariantly 
gauge invariant Weyl tensor naturally does not depend 
on the gauge choice and is always c: see Section 7 of 
PN2008. 



PN Gauge 


PN Definition 


Propagation 


PT 


PT Definition 






speed 


Gauge 




General gauge 


±PY +nU + mHU = 


c/y/n — 3 




nn + (n - 3)^x + 3(n - m)Hip = 


Chandrasekhar's gauge 


n = 3, to = arbitrary 


oo 




K + (3 - m)Hip = 


Uniform-expansion gauge 


n = 3 = to 


oo 


UEG 


K = 


Transverse-shear gauge 


n = = to 


—ic/y/3 


ZSG 


x = o 


Harmonic gauge 


n = A, m = arbitrary 


c 







Table 1. Comparison of often used gauge conditions in the PN approximation and in the perturbation theory (PT). 
The ZSG and UEG are acronyms of the zero-shear (or conformal-Newtonian, or longitudinal) gauge and the uniform- 
expansion (or uniform-Hubble) gauge, respectively; see Sections l4.1l and l4~2l The harmonic gauge in the PT is discussed 
in Section [4 .71 The gauge transformation properties and the gauge issue in the PN approach was thoroughly discussed 
in Section 6 of PN2008. The 'general gauge' condition covering most of the gauge conditions introduced in the PN 
approach was identified in Equation (210) of PN2008; here n and to are arbitrary real numbers. Notice the dependence 
of the propagation speed of the potential on the gauge choice of n; see Equation (213) of PN2008. The propagation 
speed issue (the gauge dependence of the propagation speed of potential and the real propagation speed of gravity) is 
resolved based on the Weyl tensor in Section 7 of PN2008. 



Assuming A = 
U = 0, and 



K, for the growing mode we have are 



1, 



5 = S (N) + -zU, -V • v = -6, U 



AttGqo? 
A - 



£(N) 



(31) 



where (S' N ' is the Newtonian solution with <5( N ) cx 
i 2 / 3 , t -1 ; solution for S was presented in Takada & Fu- 
tamase (1999). This implies that 2$^ + (aPi)' = 0, thus 
hi = U. The growing mode solutions valid to 1PN order 



5 \a 2 H 2 



c z 5 



(32) 



This should be compared with the Newtonian solutions 
in Equation (fTo| . Notice that only 8 has 1PN correction 
terms, whereas v and U have no PN correction to 1PN 
order. 
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3. LINEAR PERTURBATION THEORY 

We consider scalar-type perturbation in a flat 
Robertson- Walker background. Our metric convention 
is (Bardeen 1988) 



ds A = - (1 + 2a) c z dt z - 2a/3 A cdtdx l 

+a 2 [(1 + 2ip) 8ij + 27,y] dx l dx j . (33) 

We consider a flat background with vanishing pressure 
and anisotropic stress. The energy-momentum tensor is 

T ° = -{g + Sg)c 2 , T? = -a Q cu ti , Tj = 0, (34) 

where x° = cdt. To the linear order, the basic perturba- 
tion equations are (Bardeen 1988; Hwang 1991) 

k = 3Ha — 3<^ — c— 



9 A 

AttGSq + Hk + c -^cp = 0, 



A 



\2-kG 



-gav = 0, 



(35) 

(36) 

(37) 

(38) 
c c 

k + 2Hk - AnGSg + (^3H + c 2 ^ja = 0, (39) 
Sg + 3HSg - q(k- 3Ha + *v) = 0, (40) 



a" 

-X + -Hx~<p-a> = 0, 



1 . c l 

— (av) a = 0, 

a a 



(41) 



where x = a P + a2 i/ c \ K is the perturbed part of the 
trace of extrinsic curvature (minus of perturbed expan- 
sion scalar of the normal- frame vector field). These equa- 
tions are presented without taking the gauge conditions. 

Under the gauge transformation x a = x a + £, a (x e ) we 
have (Bardeen 1988; Hwang 1991) 




The gauge conditions include the hypersurface or slic- 
ing (temporal gauge) condition and congruence (spatial 
gauge) condition. Our equations are arranged using only 
spatially gauge invariant variables; for example, x is a 
spatially gauge-invariant (being independent of £ z ) com- 
bination which is the same as x m 7 — spatial gauge 
condition (Bardeen 1988). Up to this point is a review 
of the linear perturbation theory. 

Comparing Equations ([33| and (f34|) with Equations (fTJ 
and (|3]) we have 



c l c 4 



<P = ~nV, x.i = -kaPi, 
c c? 



l-flp* +2V + ZHU 



1 

p. 



(43) 



where the left and the right hand sides are perturbation 
variables in perturbation theory (PT) and in the 1PN 
approximation (PN), respectively. Notice that we have 



where 



c 

v X, 

a 



U = -c z a 



X' 



a X, 

c 



(44) 



(45) 



are gauge-invariant combinations; v x is our notation of 
gauge-invariant combination made of v and x which is 
the same as v in the zero-shear gauge % = (Hwang 
1991). To 1PN order ip, 6, Vi are gauge invariant be- 
cause we can set £( 2 ' = in the PN approximation. 
Notice that gauge transformation properties of the vari- 
ables in the PN approximation differ from the ones in 
the perturbation theory. 

In terms of the ADM (Arnowitt-Deser-Misner) and the 
covariant notations we have (Ehlers 1961; Arnowitt et al. 
1962; Ellis 1971, 1973; Bertschinger 1996) 



N = 1 — a, a, = a h R {h) = -4 



A 



7<P, 



K\ = ~ (-3H + k) , Ka - \ ViVi - ls tJ A^j X , (46) 

where TV, a c , RS h \ Kl, and Kij are the lapse function, 
acceleration vector of the normal four- vector (a c = n c; bn b 
with rii = 0) , the scalar curvature of the normal hyper- 
surface, the trace of extrinsic curvature (or expansion 
scalar of the normal four-vector 9 = n a . a with a minus 
sign), and the trace- free part of extrinsic curvature (or 
the shear-tensor o~ a b of the normal four- vector), respec- 
tively. Thus the perturbed variables a, ip, k and X can be 
interpreted as the perturbed part of the lapse function 
(or the acceleration), perturbed three-space curvature, 
perturbed part of the trace of extrinsic curvature (or the 
perturbed expansion of the normal four-vector with a 
minus sign), and the shear. Thus we naturally have the 
following names for the temporal gauge conditions: the 
synchronous gauge a = 0, the uniform-curvature gauge 
ip = 0, the uniform-expansion gauge k = 0, the zero- 
shear gauge X = 0, and the comoving gauge v = 0. 

With k defined in Equation (|35|) we can show that 
Equations (|40|) . (|4T|) and ([39j) exactly reproduce Equa- 
tions (fTTll - (fr9|l . respectively. Equation (|37|) gives 



AirGgav, = (U + HU 



(47) 



which gives Equation pT|) . 

From Equations (l39l) - (|4T|) and Equation (j4T|) . respec- 
tively, we have 



S + 2HS - AnGgS = -3 



(av)' 



2-Q 

a 



(48) 
(49) 



Notice that these forms of equation are valid without 
taking the gauge conditions. Although Equation f|48jl is 
in the same form as Eq. (|25[) with an identification a = 
—U/c 2 , we will show that only in certain hypersurface 
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conditions specified in the perturbation theory the two 
equations become identical. 

4. 1PN LIMITS OF LINEAR PERTURBATIONS 

4.1. Zero-shear hypersurface 

The zero-shear hypersurface condition takes x = 0. 
From Equation (|38]l. Equations ([35]), <(36j) and (J3H1), and 
Equation (|41l) . respectively, we have 



5 X + 2ff<5 x - 4vrG^ x = -3 



Hoe, 



a 

2-a, 

a ; 



c 2 —-'SH< 



1 . c 2 

— atv = — a v . 

a v x; a x 



3iM x = 4irGg6 



X' 



By identifying 



S = 5- 



X' 



U = 77 



2 X> 



(50) 
(51) 
(52) 

(53) 



Equation (p5"01 becomes Equation (1231) . 

We can derive closed form equations for 8 X , v x and a x . 
From Equations ([55 ]) -(|iT |) we have (Hwang & Noh 1999) 



S x + 2iJ<S x - 8ttGqS x = 



3H 2 + 6H + c 



2 A 



3H 2 



HS 



AttGq 



A 



(av x )" 
(aa x )' 



3a 2 

H (av x ) — AirGg (av x 



= o, 



i 



c 2 A \ 
\2txGqo 1 ) 



2H (aa x )' - AirGg (aa x ) = 0. 



(54) 

(55) 
(56) 



From Equation (|38p we have <p x = —a x . In the sub- 
horizon limit Equation (|54[) properly reproduces Equa- 
tion (|26p to the 1PN limit. More remarkable is the case of 
Equations l|55[) and (|56p which exactly reproduce Newto- 
nian equations in Equations (|14[) and (1151) which are valid 
to 1PN order: see Equations (1271) and (|28|) . Equations 
f5~4"]) - ((5"6"]) . are valid in fully relativistic situation and in 
the general scale. The 1PN approximation corresponds 

2 a 

to the sub-horizon limit with a c ag2 3> 1. 

The growing mode exact solutions for A = = K are 
(Hwang 1994) 



fx 



5aH 



C. 



(57) 



A complete set of exact solutions in the presence of K 
and A is presented in the Table 1 of Hwang (1994). Com- 
paring with Equations (fT6"|) and (|3^|) . in the small-scale 
limit we have 



5 X = 5^\ -Vv x =v^ 



PN) 



-fx 



>■ 



(58) 



4.2. Uniform- expansion hypersurface 

The uniform-expansion hypersurface condition takes 
k = 0. From Equation (|48|) . Equation (|39)) . and Equation 



HWANG 
(|4T|) . respectively, we have 



5 K + 2HS K - 4ttGqS k = -3 ( Ha K + 2-a K ) , (59) 



c z — +3H)a K = 4irGg6 K . 



(av K )' = —a K . 
a a 



By identifying 



(60) 
(61) 

(62) 



Equation (|59[) becomes Equation (|25[) . 

We can derive closed form equations for S K , v K and a K . 
From Equations ([35])-(|4l]) we have (Hwang & Noh 1999) 



U = --~a K , 



2H5 K - A-kGq8 k = 



3H 2 /a 
H \3H + c 2 A 



(av K ) + H (av K ) - AirGg (av K ) = —r- 

3H + c Zl % 



H + 3HHJ (av K )' + H 2 (av K ) 
(aip K ) +2H(aip K ) - AirGg {a(p K ) 
I / 3a 3 H 2 H \ 




3H + c 2 4 



,(63) 

(64) 

(65) 
(66) 



In the sub-horizon limit Equation (|63[) properly repro- 
duces Equation ([26| to the 1PN limit, whereas Equa- 
tions ((64)) and (1651) produce only the Newtonian equation. 
Equations (|63p ~ (|66p . however, are valid in fully relativis- 
tic situation and in the general scale. From Equations 
(136]) and ((39]), respectively, we have 



A 



(y3 K = -AhGq5 k , 



c z —+3H\a K = AttGq5 k , 



(67) 



either of which can be regarded as the Poisson's equation. 

The growing mode exact solutions for A = = K are 
(Hwang 1994) 



S K = 



c 2 A - 12irGga 2 
2 c 2 c 2 L 



5 aH c 2 A - 12irGga 2 
inGga 2 



5a 2 H 2 
C. 



C, 



2 A 



c 2 A 



c 2 A - 12nGga : 



:<p K . (68) 



In the large-scale limit, we have tp K = C. Comparing 
with Equations (1161) and (1321) . in the sub- horizon limit 
we have 



(69) 
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and 



-Vt)„=vW, a K = -<p K = -^U. (70) 

cr 

only to Newtonian order. 

4.3. Comoving hypersurface 

The comoving hypersurface condition takes v = 0. 
From Equation (|4ip this implies a = which is the syn- 
chronous hypersurface condition. From Equations (|39l) 
and (gU) w e have 



5„ + 2H5 V - 4nGgS v = 0, 

(a 2 K v ) + H (a 2 n v ) — inGg (a 2 n v ) 



0. 



(71) 
(72) 



Thus, in this hypersurface condition equation for 5 V is 
exactly the same as the Newtonian order equation with 
no PN correction. With an identification 



1 



-V- v, 



(73) 



equation for Kjy IS exactly the same as Newtonian one 
in Equation ([L4")) with no PN correction. Equations (|71l) 
and (|72[) are valid in the fully relativistic situation and in 
the general scale. From the gauge transformation prop- 
erties in Equation ([42) we have 






f m 




av = 




- —z \av K 

a 







Equations ([3"5), ([37)1 and (gT} give 
W = 0. 



(75) 



The coefficient of growing solution C is introduced based 
on the powerful conservation property of ip v as tp v = C. 
From the gauge transformation properties in Equation 
we have 



r . 3aH 3aH 
5 v = 5 +-#- v = 6 x + -gT v * 

3aH 



3aH 



Using this relation we can show that Equations ([50) . (|59[) 
and ((5(1)) are consistent with Eq. (|7T|) . In other words, by 
taking a combination 5 X + (3aH / c 2 )v x and using Equa- 
tions (|50l) and (|55l) we have equation for (5„ in Eq. ([7T]) . 
and similarly for S K and 5 V in Equations ([59) and ([80). 
The growing mode exact solutions for A = = AT are 



2 c 2 A 



5a 2 ff 



C, ^=C. (77) 



8 -- 2 -^C k 
5 a z H z 

Comparing with Equations (fT"6) and (|32[) . in the sub- 
horizon limit we have 



S V = S^\ - K „ = i V .v( 1PN >, 
a 



to 1PN order, and 



a v = 0, 7^ — U, 



(78) 



(79) 



even to Newtonian order. 



Notice that 5 V and k v exactly reproduce Newtonian 
density and velocity perturbations without PN correc- 
tion; this is true even to the second-order perturbations 
(Hwang & Noh 2004). Despite such remarkable identifi- 
cations, by lacking proper correspondence with PN ap- 
proximation, and lacking proper gravitational potential 
(a v — 0), it is unclear whether 5 V and be properly 

interpreted in the Newtonian context: see Section [SJ and 
for our suggestion, see below Equation ([93) . 



4.4. Uniform- curvature hypersurface 

The uniform-curvature hypersurface condition takes 
tp = 0. From Equations ||35J), ^3TJ) and (|59 > (|1T ) . and 
Equations ([31)) and ([35]), respectively, we have 



2HL - 4irGg8 v = m 2 a„ 



4nGg 
H 




(80) 
(81) 



By combining these two equations, to the 1PN order we 
have 



2H5 V - AttGqSu 



UirGga 2 
c 2 A 



H5 U 



(82) 



Compared with Equation (1261) it coincides in Newtonian 
order but differs in 1PN order. 

From Equations p5]). <j37j) and (glj we have (aHv v )' = 
and (a 3 H 2 a)- = 0, thus 



v w oc 



— i 

ag 



1 



(83) 



with no decaying mode. From Equations (I36[) and ([39 
we can derive 



ff 2 



3H 



r-2 A 



= 0. 



(84) 



The growing mode exact solutions for A = = K are 

- 2 3 



2 c 2 A 
3 — — ■ 



5a 2 ff 2 



C. 



a 



-2°- 

(85) 

Comparing with Equations ([TB]) and (1321 . in the small- 
scale limit we have 



only to Newtonian order, and 



— Vw v 7^ v 
even to Newtonian order. 



(N) 



OL v 7^ ?U, 



(86) 
(87) 



4.5. Uniform- density hypersurface 

The uniform-density hypersurface condition takes Sg = 
0. This is perfectly allowed hypersurface condition in a 
time varying background. Although the density pertur- 
bation in this hypersurface is null by definition, physical 
information like physical density perturbation is not lost 
and is hidden in other combination of variables in this 
gauge like 5 V = 6 — 3(aH/c 2 )v = —3(aH/c 2 )vs- 



NOH & HWANG 



4.6. Synchronous hypersurface 

The synchronous hypersurface condition takes a = 0. 
In the pressureless case Equation (|41[) gives v oc a -1 
which is the pure gauge mode: in Equation (I42[) a = = 
a leaves £° = ^°(x), and gives Wgaugc mode oc £°/a cx a -1 . 
By setting the gauge mode u = 0, it becomes the same 
as taking the comoving gauge. 

4.7. Harmonic hypersurface 

The harmonic or de Donder hypersurface condition 
takes perturbed part of g bc T% c = 0. The temporal com- 
ponent gives a + Ha + n = which is quite a bad choice 



leaving two remnant gauge modes: see the Appendix in 
Hwang (1993). In the PN approximation this becomes 
^P*u + 4E7 + 4HU = with the well known propagation 
speed c; however, we note that as presented in Table 1 
the propagation speed of potential depends on the PN 
gauge choice. In Table 1 we presented a combination of 
perturbation variables which reproduces this PN defini- 
tion of the harmonic gauge with a general HU term. 

5. POST-NEWTONIAN/PERTURBATION THEORY 
CORRESPONDENCE 

We summarize correspondences between the PN ap- 
proximation and the perturbation theory in Table 2. 





ZSG 


UEG 


CG 


UCG 


5 IPN 


8x 








OPN Exact 






S v [2 nd ] 




OPN Sub-horizon 


S x [2 nd ] 


S K [2 nd ] 


S v [2 nd ] 


S v [2 nd ] 


v 1PN(=0PN) 


v x 




K V [2 nd ] 




OPN Exact 


v x 




K V [2 nd ] 




OPN Sub-horizon 


[2 nd ] 


v K [2 nd ] 


Ky [2 lld ] 




U 1PN(=0PN) 


a x(= -Vx) 








OPN Exact 


a x(= -fx) 








OPN Sub-horizon 


a x (=-ip x ) [2 nd ] 


a K (= -tp H ) [2 nd ] 







Table 2. The CG and UCG are acronyms of the comoving gauge and the uniform-curvature (or flat) gauge, respec- 
tively; see Sections 14.31 and 14.41 A notation [2 nd ] means that the Newtonian correspondence is valid to second-order 
perturbation in the perturbation theory (Hwang et a!. 2012). 



As the IPN correction terms have a M order factor 

c z A 

smaller than OPN terms, those become important as the 
scale approaches the horizon. Thus the PN approxima- 
tion is applicable only in the sub-horizon scale. However, 
the Newton's theory is free from the presence of the hori- 
zon scale. This is why we have divided OPN case into the 
exact (all scales) and the sub-horizon cases in Table 2. 

Our comparison of PN approximation and perturba- 
tion theory shows that the zero-shear gauge is distin- 
guished by the fact that all perturbation variables of the 
density, velocity and gravitational potential properly cor- 
respond to IPN (thus, valid in the sub-horizon scale) 
order variables. Thus 

6 X = 5^ 1PN \ -V^,v( 1PN »=v( N >, 

a x = - Vx = ~U = -±U. (88) 

Notice that for the perturbed velocity and gravitational 
potential the Newtonian (OPN) equations remain valid 
to IPN order. 

In the uniform-expansion gauge, only the density per- 
turbation has IPN correspondence whereas the velocity 
and gravitational potential perturbations have OPN cor- 
respondences only in the sub-horizon limit. 

Newtonian correspondences of all perturbation vari- 
ables (density, velocity and gravitational potential) on 
the sub-horizon scale in the zero-shear gauge and the 
uniform expansion gauge were known to the linear order 
(Peebles 1980; Hwang & Noh 1999); the same correspon- 



dences are shown to be valid even to the second-order 
perturbations (Hwang et al. 2012). 

Newtonian correspondences of linear density perturba- 
tions in various gauge conditions in the sub-horizon scale 
were pointed out in Bardeen (1980); the same correspon- 
dences are shown to be valid even to the second-order 
perturbations (Hwang et al. 2012). Only in the comoving 
gauge (to the linear order the same as the synchronous 
gauge without the gauge mode) the density perturbation 
has an exact correspondence to the Newtonian one (Lif- 
shitz 1946; Bonnor 1957; Nariai 1969; Bardeen 1980). 

In the case of the comoving gauge the relativis- 
tic/Newtonian correspondence is quite distinguishing: 
up to the second order in perturbations the perturbed 
density and velocity show exact (thus, valid in all scales) 
correspondences with the Newtonian ones (Noh & Hwang 
2004; Hwang et al. 2012); these correspondences are pos- 
sible only in the comoving gauge. However, the comoving 
gauge implies vanishing a v , thus vanishing gravitational 
potential to the linear order. Without the proper grav- 
itational potential, despite the striking correspondences 
of density and velocity, whether proper Newtonian in- 
terpretation will be available in this gauge condition is 
unclear; for our suggestion, see below. 

For the growing modes, Equations (|32| . (15T1) and ([77]l 
show that (Takada & Futamase 1999) 



§x = Sv + ^U, 8 X = S^\ 8 V = 5<M. (89) 



POST-NEWTONIAN APPROXIMATION VS. PERTURBATION THEORY 



9 



Takada & Futamase (1999) have identified S v as the 
physical density perturbation: see below their Equation 
(A18). In perturbation theory, this solution simply fol- 
lows from the relations between the gauge-invariant com- 
binations in Equation (|76[) : 6 X = 8 V — 2>{aH / c 2 )v x . From 
Equations (|36]) - (l4T]) we can derive 

Sv = Kv, (90) 

k v + 2Hk v = 4ttGqS v , (91) 
c 2 ^a x = AhGq5 v , (92) 

These equations can be compared with Equations (|T0| - 
(|12|) in the Newtonian context. Equation (|92| is the well 
known Poisson's equations relating gauge-invariant com- 
binations based on different gauges: density perturbation 
based on the comoving gauge and gravitational poten- 
tial based on the zero-shear gauge; this was presented in 
Equations (4.3) and (4.4) of Bardeen (1980). Thus, to 
the linear order we may identify 

S v = 5, k v = — V • u, a x = — »U. (93) 
a c z 

It is likely that the proper Newtonian correspondence is 
available based on such mixed use of the gauges in con- 
structing gauge-invariant combinations. Bardeen (1980) 
in his Equations (4.5) and (4.8) has identified v x instead 
of the velocity perturbation variable; as summa- 

rized in Table 2, both v x and k v show distinguished OPN 
and 1PN correspondences, but better in having the 
exact Newtonian correspondence even to the second or- 
der in all scales. 
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